Universality aspects of the d = 3 random-bond Blume-Capel model 
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The effects of bond randomness on tlie universality aspects of tlie simple cubic lattice ferromag- 
netic Blume-Capel model are discussed. The system is studied numerically in both its first- and 
second-order phase transition regimes by a comprehensive finite-size scaling analysis. We find that 
our data for the second-order phase transition, emerging under random bonds from the second- 
order regime of the pure model, are compatible with the universality class of the 3d random Ising 
model. Furthermore, we find evidence that, the second-order transition emerging under bond ran- 
domness from the first-order regime of the pure model, belongs to a new and distinctive universality 
class. The first finding reinforces the scenario of a single universality class for the 3d Ising model 
with the three well-known types of quenched uncorrelated disorder (bond randomness, site- and 
bond-dilution). The second amounts to a strong violation of universality principle of critical phe- 
nomena. For this case of the ex-first-order 3d Blume-Capel model, we find sharp differences from 
the critical behaviors, emerging under randomness, in the cases of the ex-first-order transitions of 
the corresponding weak and strong first-order transitions in the 3d three-state and four-state Potts 
models. 
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I. INTRODUCTION 

The effect of quenched randomness on the equilibrium 
and dynamic properties of macroscopic systems is a sub- 
ject of great theoretical and practical interest. It has been 
known that quenched bond randomness may or may not 
modify the critical exponents of second-order phase tran- 
sitions, based on the Harris criterion [ll, 0- It was more 
recently established that quenched bond randomness al- 
ways affects first-order phase transitions by conversion 
to second-order phase transitions, for infinitesimal ran- 
domness in d = 2 0, H and after a threshold amount of 
randomness in d > 2 [4|, as also inferred by general ar- 
guments ^] . A physically attractive description for these 
effects has been provided by the Cardy and Jacobsen [G\ 
mapping between the random-field Ising model and the 
large q {d = 2) q-siaXe Potts model and the conjectured 
extensions of this to higher dimensions. The first nu- 
merical verification of the Cardy- Jacobsen conjecture has 
been presented recently in Ref. [7[. Furthermore, this 
rounding effect of first-order transitions has now been 
rigorously established in a unified way in low dimensions 
{d < 2) including a large variety of types of random- 
ness in classical and quantum spin systems ;»]. These 
predictions [3, Q have been confirmed by Monte Carlo 
(MC) simulations [^. Moreover, renormalization-group 
calculations on tricritical systems have revealed that not 
only first-order transitions are converted to second-order 
transitions, but the latter are controlled by a distinctive 
strong-coupling fixed point [lOj. 

Recently the present authors [11| have studied details 
of the above mentioned expectations for the effects of 
quenched bond randomness on second- and first-order 



phase transitions by considering the 2d random-bond ver- 
sion of the Blume-Capel (BC) model. Our studies have 
provided strong numerical evidence, using an efficient im- 
plementation of the Wang-Landau algorithm '12| , clari- 
fying these effects in 2d systems. In particular it was 
shown that, the second-order phase transition, emerging 
under random bonds from the second-order regime of the 
pure 2d BC model, has the same values of critical expo- 
nents as the 2d Ising universality class, with the effect 
of the bond disorder on the specific heat being well de- 
scribed by double-logarithmic corrections, our findings 
thus supporting the marginal irrelevance of quenched 
bond randomness. Furthermore, the second-order tran- 
sition, emerging under bond randomness from the first- 
order regime of the pure 2d BC model, was shown to be- 
long to a distinctive universality class with v — 1.30(6) 
and f3/i' = 0.128(5). These results pointed to the exis- 
tence of a strong violation of universality principle of crit- 
ical phenomena, since these two second-order transitions, 
with different sets of critical exponents, are between the 
same ferromagnetic and paramagnetic phases. 

In the current MC study, yielding also very accurate 
information, we wish to continue our investigations on 
the 3d BC model. We consider again a bimodal form of 
quenched bond randomness and we find different critical 
behaviors of the second-order phase transitions, emerg- 
ing from the first- and second-order regimes of the pure 
model. Hence, in the present paper we offer additional 
evidence for the anticipated strong violation of universal- 
ity. We will also constructively compare our findings to 
the existing literature on the effects of randomness in 3d 
systems, for which the pure versions undergo first- and 
second-order phase transitions. In the present MC study, 



we employ a different numerical scheme based on cluster 
algorithms and parallel tempering (PT) and, as we will 
discuss in detail, this practice is a very good alternative 
for second-order phase transitions in disordered systems. 

The relevant literature concerns mainly the 3d Ising 
model with quenched randomness, a clear case of the 
Harris criterion. The general random model has been ex- 
tensively studied using MC simulations |13l - [25j but also 
field theoretical renormalization group approaches |26l - 
[28]. The diluted model has been treated in the low- 
dilution regime by analytical perturbative renormaliza- 
tion group methods |29l - l31| and a new fixed point, in- 
dependent of the dilution, has been found. However, 
the first numerical studies suggested a continuous vari- 
ation of the critical exponents along the critical line, 
but it became clear, after the work of Ref. [T3I, that 
the concentration-dependent critical exponents found in 
MC simulations are the effective ones, characterizing the 
approach to the asymptotic regime. Nowadays, the ex- 
tensive numerical investigations of Ballesteros et al. [ij] , 
Berche et al. [2^, and Fytas and Theodorakis [2J| for 
the site-, bond-diluted, and random-bond versions of the 
model, have cleared out the doubts by presenting con- 
crete evidence that the critical behavior of the 3d Ising 
model with quenched uncorrelated disorder is controlled 
by a new random fixed point, independent of the disor- 
der strength value and the way this is implemented in 
the system. Moreover, both site- and bond-diluted Ising 
systems have been also studied by Hasenbusch et al. [20| 
in a high-statistics MC simulation, confirming universal- 
ity and providing very accurate estimates of the critical 
exponents. 

However, although the evidence of the existence of a 
unique universality class in the 3d random Ising model 
(RIM) is very strong, a full characterization of the types 
of models and the types of disorder that lead a transition 
into this class is still a challenging question. Such an in- 
teresting case is the ± J Ising model at the ferromagnetic- 
paramagnetic transition line, for which Hasenbusch et 
al. [2^ have shown that it belongs also to the RIM uni- 
versality class. Another relevant 3d random model is 
the site-diluted q — i Potts model which, in its pure ver- 
sion, is known to have a weak first-order phase transition. 
This model was studied by Ballesteros et al. [34I and 
it was clearly verified that the emerging under random- 
ness second-order transition gives a definitively different 
magnetic exponent 77. An analogous study on the 3d 
site-diluted (7 = 4 Potts model which, in its pure version, 
is known to exhibit a strong first-order phase transition 
has also produced results consistent with a further dis- 
tinctive universality class [3j] . The present study of the 
3d random-bond BC model gives us the opportunity of 
a direct comparison with the above relevant 3d cases, by 
considering the emerging second-order transitions stem- 
ming from both its Ising-like continuous and its strong 
first-order transitions, making thus our attempt more ap- 
pealing for understanding better the effect of disorder in 
3d systems. Finally, it is to be noted that, although there 



is a number of recent papers [3J-|37| dealing with the ef- 
fects of randomness on the 3d ferromagnetic BC model, 
the present study appears to be the first attempting to 
clearly distinguish between the universality aspects of the 
ex-second- and ex-first-order regimes of the 3d random 
BC model. 

The rest of the paper is laid out as follows: In the 
following Section we define the model (subsection III A[) 
and we give a detailed description of our numerical ap- 
proach (subsection [HB]) utilized to derive numerical data 
for large ensembles of realizations of the disorder dis- 
tribution and lattices with linear sizes within the range 
L S {8 — 44}. In subsection III Bl the finite-size scaling 
(FSS) scheme is described. In Sec. HIT] we discuss the 
effects of disorder on the second-order phase transition 
regime of the model. Respectively, in Sec. IIVI we illus- 
trate the conversion, under random bonds, of the origi- 
nally first-order transition of the model to second-order 
and we present full details of our FSS attempts, that 
strongly point to a new distinctive universality class. Our 
conclusions are summarized in Sec. IIVI 



II. MODEL, SIMULATION, AND FINITE-SIZE 
SCALING SCHEMES 

A. The random-bond Blume-Capel model 

The pure BC model [3^, [33| is defined by the Hamil- 
tonian 



Hp = -Jj2 s.s.+AjI^?. 



(1) 



<ij> 



where the spin variables Si take on the values —1,0, 
or -fl, < ij > indicates summation over all nearest- 
neighbor pairs of sites, and J > the ferromagnetic ex- 
change interaction. The parameter A is known as the 
crystal-field coupling and to fix the temperature scale we 
set J — 1 and fcs = 1 . 

This model is of great importance for the theory of 
phase transitions and critical phenomena and besides the 
original mean-field theory [38, 39j, has been analyzed by 
a variety of approximations and numerical approaches, 
in both d = 2 and d = 3. These include the real space 
renormalization group, MC renormalization-group cal- 
culations [40|, e-expansion renormalization grou ps l4l| . 
high- and low-temperature series calculations [42|, a 
phenomenological FSS analysis using a strip geome- 
try m, nil , and of course MC simulations [Hi l45l - [5l| . 
As mentioned already in the introduction the phase di- 
agram of the model consists of a segment of continuous 
Ising-like transitions at high temperatures and low val- 
ues of the crystal field which ends at a tricritical point, 
where it is joined with a second segment of first-order 
transitions between {At,Tt) and (Ao,T = 0). For the 
simple cubic lattice, considered in this paper, Aq = 3. 
The location of the tricritical point has been estimated 
by Deserno [4^ using microcanonical MC approach as 



the point [At,Tt] = [2.84479(30), 1.4182(55)], and to a 
better accuracy by the more recent estimation of Deng 
and Blote ^ as [At,Tt] = [2.8478(9), 1.4019(2)]. 

In the random-bond version of the BC model, stud- 
ied in this paper, the bond disorder follows the bimodal 
distribution 
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where the parameter r reflects the strength of the bond 
randomness (ratio of interaction strengths in a fixed 
50%/50% weak/strong bonds mixing). The resulting 
quenched disordered version of the model reads now as 



H 



<ij> 



JljS'iSj 'T iA ^ ^ 



si 



(3) 



Since we aim at investigating the effects of bond ran- 
domness on the segments of continuous Ising like and 
first-order transitions of the original pure model, we have 
to choose suitable values of the crystal field A and the 
disorder parameter r. For the first-order regime, we are 
looking for a convenient strong disorder combination, 
strong enough to convert the originally first-order transi- 
tion to a second-order one. The value A = 2.9 is a conve- 
nient choice, since it is well inside the region of first-order 
transitions (see the above mentioned location for the tri- 
critical point). At this value of A we have numerically 
verified that the disorder strength r = 0.5/1.5 — 1/3 
is strong enough to convert, without doubt, the original 
first-order transition to a genuine second-order transi- 
tion. Thus, wc will consider, in the sequel, the cases 
(A,r) = (1,1/3) and (A,r) = (2.9,1/3), which corre- 
spond to the ex-second- and ex-first-order regimes. These 
cases serve very well our purposes to study the univer- 
sality aspects of the emerging under bond-randomness 
second-order phase transitions of the model. 



Review^s of Monte Carlo and finite-size scaling 
schemes 



The accuracy of MC data may be decisive for a success- 
ful FSS estimation of critical properties and the proper 
selection of an algorithm is the basic requirement for the 
generation of accurate MC data. Thus, over the years, 
the numerical estimation of critical exponents has been 
a non-trivial exercise, even for the simplest models, such 
as the Ising model. An appropriate algorithm close to a 
critical point has to overcome the well-known effects of 
critical slowing down, and some otherwise excellent and 
exact algorithms, such as the Metropolis algorithm |53 |. 
which work adequately far from the critical point, be- 
come inefficient close to it. It is also well known that, 
an approach via importance sampling, close to a second- 
order phase transition, requires appropriate use of cluster 



algorithms that can efficiently overcome the critical slow- 
ing down effects. Wolff-type algorithms |53l - [55{ belong to 
this category, are easy to implement, and also very effi- 
cient close to critical points. 

Since, in the present paper, we wish to simulate the 
random-bond BC model close to critical points, the im- 
plementation of the Wolff algorithm could be a suitable 
alternative. However, for the BC model the Wolff algo- 
rithm can not be used alone, because Wolff steps act only 
on the non-zero spin values. A suggested practice is now a 
hybrid algorithm along the lines followed by Ref. [4^ . An 
elementary MC step (MCS) of this hybrid scheme con- 
sists of a number of Wolff steps (typically 5 Wolff-steps) 
followed by a Metropolis sweep of the lattice, whereas 
the usual MCS for the Metropolis algorithm is just one 
lattice sweep. The combination with the Metropolis lat- 
tice sweep is dictated by the fact that the Wolff steps act 
only on the non-zero spin values. The proposed hybrid 
algorithm is, of course, one of several alternatives for the 
study of the BC model. We note in passing that, in our 
previous analogous studies of the square lattice pure and 
random-bond BC model [llj we have been using a sophis- 
ticated two-stage Wang-Landau approach, very efficient 
for studying at the same time characteristics of the first- 
and second-order regimes of the model. 

However, since the present study concerns only the 
behavior close to the emerging, under bond random- 
ness, critical points, the described hybrid approach, com- 
bined with a PT protocol, which will be detailed be- 
low, is much more appealing and its convergence may 
be easily checked. Thus, we have simulated the random- 
bond BC model on the simple cubic lattice for the cases 
(A,r) = (1,1/3) and (A, r) = (2.9,1/3) and for lattice 
sizes in the range L = 8 — 44, by implementing the hy- 
brid approach described above, suitably adapted to the 
present disordered model. As is always the case, the sim- 
ulation task for a disordered system, such as the present 
random-bond BC model, is by far more demanding than 
its pure counterpart, since one has to sum over disorder 
realizations. Furthermore, for the random-bond model 
the simple Wolff algorithm cannot be applied and a suit- 
able generalization is necessary. This generalization is 
a straightforward adaption of the Wolff algorithm by 
introducing two different probabilities for putting links 
between sites; corresponding to the weak and strong 
nearest-neighbor ferromagnetic interactions. It is then 
rather easy to show that, the generalized cluster algo- 
rithm has all the properties of the original Wolff algo- 
rithm. 

In order to extract information on the critical prop- 
erties of the random-bond BC model, we are going to 
apply the ideas of FSS. For the estimation of the critical 
temperatures and the corresponding critical exponents, 
one has to generate MC data to cover several finite-size 
anomalies of the finite systems of linear size L. The above 
described hybrid approach has to be carried over to a 
certain temperature range depending on the lattice size. 
These temperatures may be selected independently or it 
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FIG. 1: (color online) Behavior of moving averages for the 
energy per site {E/N}^, in panel (a), and the order param- 
eter {M)^, in panel (b), at a selected temperature, close to 
the critical. The behavior of the four simulation schemes is 
illustrated and the definition of elementary steps is discussed 
in the text. 



may also be convenient, and accurate, to implement a PT 
protocol. Such a PT protocol, based (mainly) on tem- 
perature sequences corresponding to an exchange rate 
0.5, was tested for both cases considered in this work, 
i.e. (A,r) = (1,1/3) and (A, r) = (2.9,1/3), and was 
found to be very accurate. We should point out here that, 
the proposed PT approach is very close to the practice 
suggested recently in Ref. |56|. Furthermore, a similar 
approach has been followed recently for a very success- 
ful estimation of critical properties of the Ising and the 
S = 1 model (BC) on some Archimedean lattices [57| . 
Appropriate temperature sequences, for the application 
of a hybrid PT protocol were generated via short pre- 
liminary runs. Dur ing these preliminary runs a simple 
histogram method [53|, [5^ was applied to yield energy 
probability density functions and from these, the appro- 
priate sequences of temperatures were easily determined, 
following the recipe of Ref. [56j. The preliminary runs 
were applied to only one disorder realization, but is was 
checked that the resulting temperature sequences were 
producing, to an excellent approximation, a constant ex- 
change rate of 0.5 for the PT protocol, when applied to 
any other different disorder realization. This is of course 
an intrinsic, and very convenient property of the present 
random-bond model. 

The above MC scheme was carefully tested almost for 
all lattice sizes before its implementation for the gen- 
eration of MC data necessary for applying FSS. These 



tests included two basic ingredients: the first was the 
estimation of the MC times necessary for equilibration 
and thermal averaging process applied to a particular 
disorder realization and the second involved the observa- 
tion of running sample-averages of some thermodynamic 
quantity, such as for instance, the magnetic susceptibil- 
ity, at some particular temperature close to the critical 
one. These tests, not only provided convincing evidence 
for a successful determination of equilibrium results, but 
also highlighted the way towards an efficient and opti- 
mum scheme. For instance it was checked that for a 
rather small lattice size, such as L — 16, the sample 
averages over 150 disorder realizations of all thermody- 
namic quantities were identical, irrespective of the al- 
gorithm used (simple Metropolis or hybrid algorithm), 
provided one uses reasonable MC times for equilibration 
and thermal averaging. This could be explained as a 
result of effective cancellation of statistical errors from 
the sample-averaging process. Noteworthy, that sample- 
to-sample fluctuations are much larger than any usual 
statistical errors in a proper MC thermal averaging pro- 
cess. Of course, for a particular disorder realization the 
difference in the outcome of different algorithms may be 
truly spectacular. 

In fact, we can observe the superiority of the hybrid 
approach, over a simple Metropolis scheme [52|, in Fig.[T] 
This figure is constructed by using moving averages for 
the energy per site {{E/N)t) and the order parameter 
{{M)t) for a lattice size L = 40 at a selected tem- 
perature, close to the critical temperature, for the case 
(A,r) = (1, 1/3) of the random-bond BC model. As it 
can be seen from this illustration, the Metropolis algo- 
rithm suffers from very strong fluctuations and so does 
the corresponding PT protocol (applied in a tempera- 
ture sequence including the temperature shown) using 
this algorithm. They both follow a very slow approach 
to equilibrium and only with the help of extensive sam- 
pling one could expect results of reasonable accuracy. On 
the other hand, the hybrid approach converges very fast 
to equilibrium and produces quite accurate results. The 
corresponding PT hybrid algorithm appears to further 
improve this good behavior, which is a reasonable expec- 
tation, since the PT approach should be in general more 
effective for disordered systems. Note also that. Fig. [T] is 
constructed by using one particular disorder realization 
and the dashed straight lines in the panel of the figure 
indicate the final average values of four independent PT 
hybrid runs of the same realization, whereas the fluctuat- 
ing lines correspond to one single run for each algorithm. 
A careful examination of this illustration shows that even 
a single run for the thermal averaging process is very ac- 
curate for the hybrid and the PT hybrid algorithms and 
for the presented case, a MC time of the order of 3 x A^, 
is adequate and optimum for the thermal process. The 
corresponding equilibration MC times are always shorter 
than this. There is no doubt that, the use of the hybrid 
algorithm is essential for the study of the present model. 

Let us review now details from the FSS tools used 



throughout the paper, for the estimation of critical prop- 
erties of the disordered system. The outhne below follows 
our previous studies on random-bond models and is very 
close to the traditional standard tools of the theory of 
FSS. When studying a disordered system, a large num- 
ber of disorder realizations has to be used in the sum- 
mations in order to obtain good sample-averages of any 
basic thermodynamic quantity Z, which is a usual ther- 
mal average of a single disorder realization. The sample 
(or disorder) averages may be then denoted as \Z\av and 
their finite-size anomalies respectively as [Z]^^. These 
(disorder-averaged) finite-size anomalies will be used in 
our FSS attempts, following a quite common practice [33[ 
and their temperature locations will be denoted by T[z]'^ ■ 
Thus, we also follow the conventional numerical approach 
that appears to work well for the present model, at both 
the ex-second- and ex- first-order regimes, at the disorder 
strength values considered. However, in an alternative 
approach 18] one may consider individual sample depen- 
dent maxima (anomalies) and the corresponding sample 
dependent pseudocritical temperatures. This alterative 
route is far more demanding computationally, but the 
corresponding FSS analysis may be more precise, and ad- 
ditional useful information concerning the properties of 
disorder averages becomes available. Our study addresses 
only exponents describing the disorder-averaged behav- 
ior and we do not use a FSS analysis based on sample 
dependent pseudocritical temperatures. For disordered 
systems one can make a clear distinction between typical 
and averaged exponents [58, 59]. An extreme example of 
this case with very pronounced differences in the corre- 
sponding V exponents has been provided, and critically 
discussed, by Fisher [531 O"^ the random transverse-field 
Ising chain model. Finally, disordered critical phenom- 
ena are known to display in general multicritical expo- 
nents [60J, |6l| and we shall return to this interesting point 
in our last section, when discussing the violation of uni- 
versality in the present model. 

The number of disorder realizations and the selection 
of temperatures may influence the accuracy and suitabil- 
ity of the MC data from which the locations of the finite- 
size anomalies are determined, by fitting a suitable curve 
(for instance a fourth-order polynomial) in the neighbor- 
hood of the corresponding peak. Since we are imple- 
menting a PT hybrid approach, based on temperatures 
corresponding to an exchange rate 0.5, we are selecting 
certain temperature sequences consisting of a number of 
(say 3 or 5) different temperatures, averaging then over 
a relatively moderate number (~ 100) of disorder real- 
izations. This yields a number of (3 or 5) points of the 
averaged curves [Z]av This set of points may not be 
adequately dense and will not cover the ranges of the 
peaks of all thermodynamic quantities of interest. There- 
fore, this procedure is repeated several times (depending 
on the linear size L) by using new sets of temperatures 
which are translated with respect to the previous sets. 
These translations should be carefully chosen so that a 
final dense set of points, suitable for all finite-size anoma- 



lies, is obtained. This corresponds, on average, to a very 
large number of realizations, since for each PT hybrid 
run, a different set of disorder realizations was used. We 
have found this practice very convenient, efficient, and 
most importantly quite accurate. Thus, we were able to 
describe the averaged finite-size anomalies of the system 
with high accuracy. 

In order to estimate the critical temperature, we follow 
the practice of simultaneous fitting approach of several 
pseudocritical temperatures [Ilj. From the MC data, 
several pseudocritical temperatures are estimated, corre- 
sponding to finite-size anomalies, and then a simultane- 
ous fitting is attempted to the expected power-law shift 
behavior T[z]*^ =Tc + bz- L'^^". The traditionally used 
specific heat and magnetic susceptibility peaks, as well 
as the peaks corresponding to the following logarithmic 
derivatives of the powers n — 1, 2, and rt = 4 of the 
order parameter with respect to the inverse temperature 
K=l/T [13, 
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and the peak corresponding to the absolute order- 
parameter derivative 
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will be implemented for a simultaneous fitting attempt of 
the corresponding pseudocritical temperatures. Further- 
more, the behavior of the crossing temperatures of the 
fourth-order Binder's cumulant [63], and their asymp- 
totic trend, will be observed and utilized for a safe esti- 
mation of the critical temperatures. 

The above described simultaneous fitting approach 
provides also an estimate of the correlation length ex- 
ponent ;/. An alternative estimation of this exponent is 
obtained from the behavior of the maxima of the loga- 
rithmic derivatives of the powers n = 1, 2, and n = 4 
of the order parameter with respect to the inverse tem- 
perature, since these scale as '^ L^^'^ with the system 
size [63] ■ Once the exponent ly is well estimated, the 
behavior of the values of the peaks corresponding to 
the absolute order-parameter derivative, which scale as 
^ L^'^^l^')/'^ with the system size [62], gives one route for 
the estimation of the magnetic exponent ratio /3/iy. Addi- 
tionally, knowing the exact critical temperature, or very 
good estimates of it, we can utilize the behavior of the 
order parameter at the critical temperature for the tra- 
ditionally effective estimation of the exponent ratio /3/j/ 
(Afc = M{T = Tc) - L-P/"). Summarizing, our FSS 
approach utilizes, besides the traditionally used specific 
heat and magnetic susceptibility maxima, the above four 
additional finite-size anomalies for an accurate estima- 
tion of critical temperatures and relevant exponents. 



III. EX-SECOND-ORDER REGIME: 

UNIVERSALITY CLASS OF THE RANDOM 

ISING MODEL 



We study in this Section the 3d random-bond BC 
model at the second-order regime of the phase diagram. 
We have considered the case (A,r) = (1, 1/3) and simu- 
lated lattices with linear sizes in the range L = 8 — 44. 
Below we discuss in detail our efforts to find a compre- 
hensive FSS scheme to fit the numerical data. However, 
before illustrating details of our study, it is useful to 
recall again that, according to general universality ar- 
guments, the pure BC model at A = 1 is expected to 
belong to the 3d Ising universality class and thus, the 
case studied here, should be contrasted to the relevant 
literature on the 3d RIM. From the relevant literature, 
it appears that a consensus have been achieved today 
with regard to the existence of a single universality class 
for the general 3d RIM [13, IH, Hi, [H. However, it 
is also true that the study of Berche et al. [22| on the 
3d bond-diluted Ising model, has emphasized the strong 
influence of crossover phenomena and pointed out that 
the identification of this universality class is a very diffi- 
cult task. For the case of magnetic bond concentration 
p — 0.7, studied in that paper, these authors found an 
effective exponent {l/v)f.ff = 1.52(2). Correspondingly, 
for the case p = 0.55 the effective value was found to 
be {l/i')eff — 1.46(2). This noticeable variation in the 
effective exponents indicated possible confiucnt correc- 
tions and/or crossover terms ^22]. They finally agreed 
that, the case p — 0.55 appears to be the case of least 
scaling corrections and thus, in this way, agreement with 
the value given by Ballesteros et al. [17] was established. 

The problem of slowing decaying scaling corrections 
has been discussed in detail for both these randomly site- 
and bond-diluted systems by Hasenbusch et al. [20] . and 
for the latter, the value p — 0.54(2) has been proposed 
as the value at which the leading scaling corrections van- 
ish. The above observations will be very useful later in 
this Section, when we discuss our problems in extract- 
ing a reliable estimate of this exponent, since we have 
faced similar problems. Noteworthy that, in a system- 
atic study of the ± J Ising model at the ferroniagnetic- 
paramagnetic transition line, Hasenbusch et al. [25{ have 
analyzed the effects of leading and next-to-leading scal- 
ing corrections and proposed a value of the concentration 
of the ferromagnetic lines of this model at which correc- 
tions to scaling almost vanish. The final estimate for the 
correlation length exponent determined by this study is 
almost identical with the value given by Ballesteros et 
al. [13. 

We start the FSS analysis by attempting the estima- 
tion of the exponent ratio that characterizes the diver- 
gence of the susceptibility. We assume that the (sample- 
averaged) susceptibility obeys a simple power law of the 
form [xt]q„ = b ■ L^/" . The fitting attempts are found to 
be quite stable with respect to the size range chosen. We 
attempted the ranges L = 8 - 44, 12-44, . . . , 28-44 



■ [x,]„ = 0.0633(4) L' 




15 20 25 30 

L 



35 



40 



45 



FIG. 2: (color online) Illustration of the divergence of the sus- 
ceptibility maxima by a simultaneous global fitting attempt 
on the data corresponding to several lattice-size ranges shown 
in the panel, indicated by the different line colors. 



and observed the behavior of the resulting coefficients h 
and exponents ^ /v. The estimates of the exponent vary 
slowly in the range ^ /v = 1.962 — 1.967, as we increase 
Lmin from L,-ain = 8 to L,„i„ = 28. A kind of mean 
value estimate can be obtained by a global fitting at- 
tempt simultaneously applied to the above ranges. Such 
a simultaneous attempt using the ranges L = 16 — 44, 
20 - 44, 24 - 44, and L = 28 - 44 is shown in Fig. [5] giv- 
ing an estimate ^/v = 1.963(3). The illustrated estimate 
is in full agreement with the values given by Ballesteros 
et al. PJI] and Berche et al. [22|- Finally, observing the 
over all trend of such simultaneous fitting attempts we 
shall propose as our final estimate ^/v = 1.964(4), in 
coincidence with the estimate ^/v — 1.965(10) given by 
Berche et al. [24I for the case p — 0.7 mentioned above. 

A similar simple power-law behavior was found for 
the peaks corresponding to the absolute order-parameter 
derivative which scale as ~ Ij0--P)/'^ _ Since the exam- 
ined behavior, in this case, is very similar to the above 
behavior, we shall not discuss further our efforts and 
give only our final estimate for the relevant exponent 
(1 — fi)/v = 1.022(5). Combining the above estimates 
and assuming at this point hyperscaling - in the form: 
2/v = 2(1 - 13) /v + {d- -f/v) - we find l/u = 1.54(1). In 
the following we will see that, this value is not far from 
the effective estimates determined below from the shift 
behavior of the system. 

Let us now attempt the estimation of the correlation 
length exponent via the scaling behavior of the logarith- 
mic derivatives of the powers n = 1, 2, and n = 4 of 
the order parameter with respect to the inverse temper- 
ature [see Eq. (g])]. Their behavior was observed to be 
quite successfully fitted to a stable simple power law with 
rather small variation of the effective critical exponent. 
All the estimates were close to the above value \/v = 1.54 
and for illustrative reasons we have presented one such 




FIG. 3: FSS behavior of the peaks of the logarithmic deriva- 
tives of the powers n = 1,2, and n = 4 of the order parameter 
with respect to the inverse temperature. The estimate for the 
exponent 1/v is given in the panel by applying a simultane- 
ous fitting attempt to a simple power law in the size range 
L = 16 - 44. 
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FIG. 4: (color online) Illustration of the asymptotic trend of 
the fourth-order Binder's cumulant of the order parameter for 
various temperatures close to the critical temperature. Note 
the stability of the estimates at Tc — 2.88035. 



simultaneous fitting attempt in Fig. [31 Therefore, the 
scheme appears to be in good agreement with hyperscal- 
ing. We turn now to the estimation of this exponent from 
the general shift behavior of the system. 

As is well known, the pseudocritical temperatures, cor- 
responding to several finite-size anomalies, provide a tra- 
ditional route for the estimation of the critical temper- 
ature and the correlation length exponent. A simulta- 
neous fitting attempt to a power-law shift behavior of 
the form T^z]* = Tc + bz ■ L^^/^ is often the attempted 
practice. For the present case the following tempera- 
tures were calculated. Temperatures of the peaks of the 
specific heat, magnetic susceptibility, inverse tempera- 
ture derivative of the absolute order parameter, and in- 
verse temperature logarithmic derivatives of the n — 1, 
2, and n = 4 powers of the order parameter. The data 
were fitted in the ranges L = 8 - 44 to L = 28 - 44 
and the behavior of the estimates was observed. With 
the exception of the fitting attempt corresponding to 
L = 8 — 44, which gave an estimate 1/v = 1.465(1), 
all other attempts gave fluctuating estimates in the range 
l/i/ — 1.52 — 1.54, which seems to agree well with the pre- 
vious finding. Yet, at the same time, we observed a no- 
ticeable shift of the estimated critical temperature from 
values Tc = 2.8798 to Tc = 2.88015 as we varied Lmin 
from Lmin = 8 to Lmin = 28. Here, a careful examina- 
tion of the behavior of the fourth-order Binder's cumu- 
lant of the order parameter [VT]av = [l-{M'^)/3{M^)^]av 
suggested that a larger critical temperature in the range 
Tc = 2.88015-2.88045 may be a very strong option. Fig- 
ure m is a very clear illustration of this statement, since 
it demonstrates that, at the temperature Tc — 2.88035, 
the behavior of the cumulant is almost independent of L. 

In order to better understand the shift-behavior of 



the system and present a global illustration of it, we 
fix the critical temperature to several values and calcu- 
late effective exponents. We studied their size depen- 
dence by varying Lmin as usual, from Lmin = 12 to 
Lmin = 28, and by further applying simultaneous fit- 
ting attempts using the form T^z]^^ = Tc + bz ■ L~^^^ 
with fixed critical temperatures from the set Tc = 
{2.88015, 2.88025, 2.88035, 2.88045, 2.8807}. The re- 
sulting global behavior is now illustrated in Fig. [SJ This 
figure clearly indicates that effective exponents are very 
sensitive to small temperature changes. Furthermore, 
by applying a linear fitting in the effective estimates for 
L > 16 for the case Tc = 2.88035 we obtain l/v = 1.504 
and this is indicated by the bold straight line in the panel, 
as our central estimate for the exponent. Repeating the 
same procedure for the case Tc = 2.88025 we obtain 
\/v = 1.539, whereas for the case Tc = 2.88045 we ob- 
tain the value l/iy = 1.466. We may use these remote 
estimates as upper and lower bounds respectively of our 
exponent estimation and these are indicated in the panel 
by the doted lines. For the sake of comparison we show 
also in the panel the accepted limits for the estimation 
of the critical exponent 1/v of the pure 3d Ising model 
1/ = 0.6304(13) and the 3d RIM i/ = 0.6837(53) [13 ■ 
Based on the above observations, one could anticipate 
that any estimate in the range l/i^ — 1.466 — 1.539, indi- 
cated by the dotted lines in the panel, should be accept- 
able. We point out here that, the exponent for the pure 
3d Ising model, j/ = 0.6304(13) [H] illustrated in Fig. [3 
is a moderate estimate with rather large error bounds. 
A most recent and accurate estimate is j^ = 0.63002(10), 
obtained by Hasenbusch [6^ in excellent a gree ment with 
other recent studies of the 3d Ising model [6a, [63| ■ 

Overall, the present case shares many of the problems 
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FIG. 5; (color online) A global illustration of the estimates 
of the effective exponent (l/i/)e//. The solid line drawn in 
the panel close to the value l/v = 1.586, together with the 
dashed lines, illustrate the location of the critical exponent 
range for the pure 3d Ising model. The analogous range for 
the 3d RIM is drawn close to the value 1/v = 1.463. The 
range l/f — 1.466 — 1.539 (dotted lines in the panel) around 
the heavy solid line at the value l/u = 1.504 demonstrate the 
possible asymptotic evolution of the critical exponent for the 
present model. 



encountered by Berche et al. [22| in their study of the 
3d bond-diluted Ising model for the case of magnetic 
bond concentration p = 0.7. However, the above illus- 
trations show that, despite the already mentioned prob- 
lems, our data are compatible with the general expec- 
tations of the general 3d RIM. We close this Section 
with two additional comments. Firstly, we also tried 
to observe the shift behavior by fixing the shift expo- 
nent to the value l/f = 1.463 proposed by Ballesteros 
et al. [13] ■ We then found that the corresponding si- 
multaneous fitting attempts gave values for the criti- 
cal temperatures which are slowly approaching the value 
Tc = 2.88035(10) suggesting once again that this may 
be the asymptotic critical temperature. Secondly, we ob- 
served the FSS behavior of the order parameter by fix- 
ing the critical temperature to the values from the set 
Tc = {2.88015, 2.88025, 2.88035, 2.88045}. The esti- 
mated effective exponents iS/iy, had in all cases values in 
the ranges /3/i/ = 0.51(1) - 0.52(1). These values, to- 
gether with the earlier estimate 7/1^ — 1.964(4), satisfy 
quite well hyperscaling. 



IV. EX-FIRST-ORDER REGIME: 

DISORDER-INDUCED SECOND-ORDER PHASE 

TRANSITION 

In this Section we investigate the effects of bond ran- 
domness on the segment of the first-order transitions of 
the original pure model. At the value A — 2.9, which 



FIG. 6: (color online) Illustration of the continuous behav- 
ior of the order parameter for the random-bond BC model 
at (A,r) = (2.9,1/3). Three different lattice sizes L = 12, 
16, and L = 36 are shown and the illustrated behavior is ap- 
proximately centered at the temperatures corresponding to 
the peaks of the magnetic susceptibility, indicated by the as- 
terisks and the drop lines. 



is well inside the first-order regime of the pure model, 
we have numerically verified that the disorder strength 
r = 1/3 is strong enough to convert, without doubt, the 
original first-order transition to a genuine second-order 
one. The illustration in Fig. |6] consists very clear evidence 
of this statement. In this figure we show the behavior of 
the order parameter close to the pseudocritical region for 
three different lattice sizes L = 12, 16, and L — 36. The 
three different asterisks and the respective drop lines in- 
dicate the points on these curves corresponding to the 
magnetic susceptibility peaks. These order-parameter 
curves, at the corresponding pseudocritical regions, are 
smooth and no sign of discontinuity appears as the lat- 
tice size increases to L — 36. Since the order parame- 
ter is continuous at the combination (A, r) = (2.9, 1/3), 
the disorder-induced transition is a second-order phase 
transition, which we will take as representative of the 
ex-first-order universality class. 

For the case (A, r) = (2.9, 1/3), we start the FSS anal- 
ysis by attempting the estimation of the exponent ra- 
tio that characterizes the divergence of the susceptibility. 
First, we assume that the (sample-averaged) susceptibil- 
ity obeys a simple power law of the form [xt]qi, = b-L''''^ . 
However, the fitting attempts are found to be unstable 
with respect to the size range (we attempted the ranges 
L — 8 — 44, 12 — 44, . . .) and the resulting coefficients 
b and exponents j/v are varying in a competitive way, 
indicating the need of correction terms. Therefore, we 
tried to find a resolution of this problem by introduc- 
ing correction terms and we found that a constant back- 
ground term, was actually the most effective in eliminat- 
ing the observed instability. The behavior of the corre- 
sponding fitting attempts, assuming the scaling relation 



100 


- 


1 1 


' 1 ' 

^=1.10(15) 


+ 


' 1 
087(2) 


' 1 ' 1 


1 ' 1 

/ 


- 


80 


- 












/ 


- 


S 60 


_ 










X 




_ 


H 
X 












0" 




■ 


■— 40 








rt^ 


X 


□ L>8 


• L>12 


- 


20 






^-^^ 






A L> 16 
♦ L>24 


L>20 
L>28 







- 


1 1 


1 1 




■ 1 


A L>32 

1 1 1 1 


1 i 1 


- 



20 



25 
L 



•^ 10 




• L>28: [xj,^,- 1.1 =0.0871(4) L 



25 30 35 40 45 



FIG. 7: (color online) The divergence of the susceptibil- 
ity maxima is well described by a scaling law of the form 
[xtIou = a + 6 ■ L"^' " . Illustration of a simultaneous global fit- 
ting attempt on the data corresponding to several size ranges 
shown in the panel, indicated by the different line colors. 
Note that, the simultaneous attempt gives an estimate of 
a = 1.10(15) which is in agreement, within errors, with all 
separate attempts and effectively coincides with the estimates 
of the ranges L = 8 - 44 and L = 12 — 44. 
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b ■ L'^'" , is very good and produces very 



slowly varying fitting parameters. A completely stable 
behavior was obtained by fixing the value of the back- 
ground term to be of the order ofa= 1.0 — 1.2. The 
values of tlie background term for the first two ranges 
L = 8 — 44 and L = 12 — 44, and the value obtained by 
a simultaneous global fitting attempt on tlie data corre- 
sponding to several size ranges, as shown in Fig. [71 are 
effectively the same. As indicated in the panel of this fig- 
ure the global attempt gives an estimate a = 1.10(15). In 
the same panel we give the estimates of 5 = 0.087(2) and 
the exponent j/i/ = 1.863(6). Using this value (a = 1.1) 
we illustrate in Fig.[5]the stability of the resulting scaling 
scheme by presenting, in a double logarithmic scale, the 
behavior of the two most remote ranges L = 8 — 44 and 
L = 28 — 44. The coincidence of the estimates in this 
panel appears as a guaranty of the quite accurate esti- 
mation of the magnetic exponent ratio 7/1/ ~ 1.864(12). 

However, it may be crucial for future studies to give 
here a more detailed discussion on the various corrections 
tested before adopting the above scenario. Our fitting 
attempts followed a quite common practice restricting 
the search to an expression including only one correc- 
tion term, e.g., [xrllv = b ■ 1^/" + b' ■ L^/"-". This 
restriction was unavoidable, since the fits were already 
notoriously unstable, as we increased the minimum lat- 
tice size, with the above four-parameter expression. The 
stability and quality of the fittings were then observed 
by fixing the correction exponent e to various relevant 
values and treating only the exponent 7/1^ as a free one. 
The quality of the fittings was characterized by their cu- 



FIG. 8: (color online) This figure is complementary to Fig. [7] 
and further illustrates the stability of the scaling scheme by 
presenting, in a log-log scale, the behavior of the two most 
remote ranges L = 8 — 44 and L — 28 — 44. Note that, we have 
here subtracted from the susceptibility data the estimated 
value of the background term (a = 1.1). 



mulated square deviation, x^, and the stability of the 
values of the exponent 7/1/ and the corresponding ampli- 
tudes were observed. The best sequence of fittings was 
obtained for e = "f/v, corresponding to the regular back- 
ground behavior adopted in this paper and illustrated 
in Figs. [7] and [5] This choice is unique, in the sense 
that is for sure the simplest one. Besides, all the other 
tested values of e, in the neighborhood of the leading 
and next-to-leading corrections to scaling, correspond- 
ing to the 3d RIM e = w = 0.33(3), e = 2w = 0.66(6) 
and e = UJ2 = 0.82(8) [20|], did not produce a stable 
sequence of fittings, but gave a rather pathological vari- 
ation of the estimates of the exponent j/v and the rel- 
evant amplitudes. Yet, one may also observe the stabil- 
ity and quality of the fittings by fixing both the expo- 
nent "f/v and the exponent e. In these final attempts, 
we did found cases of comparable quality and stability 
that are completely compatible with the 3d RIM univer- 
sality class. One such case is described by the expres- 
sion [xt]L = 0.0494(4) • Li-964 + 0.013(1) • Li-964-o.66^ 
and another one is described by the expression [xtIqu = 
0.052(2) ■ L^-^^'^ + 0.18(2) ■ 2,1-964-0.82^ rpj^g ^^^^^ ^^^ 

expressions and the behavior in Fig. [5] produce very close 
values (almost identical within statistical errors) in the 
range L = 8 — 44 studied in this paper. This obser- 
vation should serve also as a warning of the difficulties 
and the pathology of the fitting attempts, since the exis- 
tence of stable forms, with quite different correction ex- 
ponents, means also that a completely reliable estimation 
of corrections-to-scaling exponents may not been feasible 
even at larger lattice sizes. 

Let us now investigate the FSS behavior of the peaks 
corresponding to the absolute order-parameter deriva- 
tive which, as mentioned earlier, are expected to scale 
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FIG. 9: (color online) Illustration of a simultaneous global fit- 
ting attempt on the data corresponding to several size ranges 
shown in the panel for the maxima of the sample-averaged 
absolute order-parameter derivative. Note that in this case, 
the fitting attempts are insensitive to the size ranges used. 



FIG. 10: FSS behavior of the pseudocritical temperatures de- 
fined in the text. Estimates for the critical temperatures and 
the shift exponent Xjv are given in the panel. The stability 
of the fitting scheme is discussed in the relevant text. 



as ~ L'^^ ^^I" with the lattice size. Here, we find that, 
the corresponding maxima obey very well a simple power 
law [a(|M|)/9ii:];„ = h ■ L(i-'3)/'' without any correction 
terms. The corresponding fitting attempts are very sta- 
ble with respect to the lattice-size range, as moving from 
L = 8-44 to L = 28-44. Figure IHliUustrates in the main 
panel a simultaneous global fitting attempt on the data 
corresponding to the size ranges shown. The inset of this 
figure presents a simple fitting for the complete lattice 
range L — % — 44. The estimates for the exponent are 
almost insensitive to the used lattice-size ranges giving 
for the simultaneous global fitting (1 — /3)/i' = 0.870(5) 
and for the simple fitting L = 8 — 44, in the inset, 
(1 — /3)/i^ = 0.875(6). These results indicate the con- 
sistency of the FSS scheme and the accuracy of the nu- 
merical data. We may take as a final (quite confident 
and moderate in its error bounds) estimate, the value 
(l-/3)/i^ = 0.87(1). 

As mentioned earlier, the pseudocritical temperatures, 
corresponding to several finite-size anomalies, provide 
a route for the estimation of the critical temperature 
and the correlation length exponent. A simultaneous fit- 
ting attempt to a power-law shift behavior of the form 
T[z]*^ = Tc+bz-L'^/" is the generally suggested practice. 
Figure [TU] illustrates the shift behavior of such several 
pseudocritical temperatures. These temperatures corre- 
spond to the peaks of the following six (sample-averaged) 
quantities: specific heat, magnetic susceptibility, inverse 
temperature derivative of the absolute order parameter, 
and inverse temperature logarithmic derivatives of the 
n — \^ n — 2^ and n — A powers of the order parameter. 
The data illustrated are fitted in the range L = 12 — 44, 
which corresponds to the best fitting attempt of all tried. 
The resulting estimates of the critical temperatures and 
the shift exponent 1/v are given in the panel. However, 



in this case also, the fitting attempts for the estimation of 
the critical exponent 1/v from the data of the pseudocrit- 
ical temperatures are not completely stable. Therefore, 
some further comments and analysis are necessary. De- 
spite that, the estimation of the critical temperature is 
quite stable, with values ranging from Tc = 1.6833(2) to 
Tc = 1.6835(2), as we vary the lattice-size ranges from 
L = 8 — 44 to L = 24 — 44. However, the corresponding 
exponent estimates 1/v have a noticeable variation with 
estimates from 1.51(2) to 1.42(2). The statistical errors 
influence here the quality but also the stability of the 
fitting attempts. 

To better understand the shift-behavior we tried the 
following two assumptions. First, we fixed the critical 
temperature to be Tc = 1.6835, a value indicated also 
by a very careful examination of the behavior of the 
fourth-order Binder's cumulant of the order parameter 
(not shown here for brevity). Again, we found a similar 
variation as above with the lattice-size ranges used. The 
best fittings were obtained for the ranges L = 12 — 44 
and L — 16 — 44, giving also very close estimates for 
the shift exponent 1/v [1/v = 1.45(2)]. Subsequently, we 
tried to observe the behavior of the estimates of the criti- 
cal temperature by fixing the shift exponent to the value 
1/v = 1.438. This is the value obtained by satisfying hy- 
perscaling, given the previous estimates for "f/v = 1.864 
and (1 — P)/v = 0.87. The estimates for the critical 
temperature are all very close to Tc = 1.68345 and the 
best fitting attempt, giving also an estimate with the 
smallest error, is Tc = 1.6835, corresponding to the range 
L = 12 — 44. We have therefore accepted as our best es- 
timation of the shift behavior the values Tc = 1.6835(2) 
and l/v = 1.45(2). 

An independent estimation of the correlation length 
exponent can be obtained via the scaling behavior of the 
logarithmic derivatives of the powers n = 1, 2, and n = 4 
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FIG. 11: FSS behavior of the peaks of the logarithmic deriva- 
tives of the powers n = 1,2, and n = 4 of the order parameter 
with respect to the inverse temperature. The estimate for the 
exponent 1/v is given in the panel by applying a simultane- 
ous fitting attempt to a simple power law in the size range 
L = 12- AA. 



of the order parameter with respect to the inverse tem- 
perature [see Eq. @]. Their behavior was observed to 
be quite stable and consistent with the above estimation 
from the shift behavior. We tried here to vary again the 
ranges from L = 8 — 44 to L = 28 — 44. Depending 
on Lmin the estimated effective exponents vary in the 
range \/v = 1.445 — 1.415. In Fig. [Tl]we illustrate such 
an estimation by a simultaneous fitting attempt in the 
range L — 12 — 44. As it can be seen from this figure 
the estimate is 1/v = 1.440(2). Therefore, from this FSS 
scheme one should conclude that 1/v = 1.430(10) and 
combining with the above shift behavior, we should re- 
gard 1/v = 1.440(10) a very decent proposal that is now 
in full agreement with the exponent value 1/v = 1.438 
obtained by satisfying hyperscaling, given the previous 
estimates for j/v = 1.864 and (1 - l3)/v = 0.87. 

Finally we give an outline on the behavior of the or- 
der parameter at the estimated critical temperature. We 
computed the finite-size values of the order parameter 
at the temperature Tc = 1.6835 from the corresponding 
(sample-averaged) order-parameter curves. In Fig. [12] we 
apply a simple power-law estimation for the exponent ra- 
tio [3/v, using again the size range L = 12 — 44. This es- 
timation gives a critical exponent ratio /S/v = 0.566(5) in 
excellent agreement with the value /3/v = 0.568, which is 
obtained by satisfying hyperscaling, given the estimates 
-f/v = 1.864 and (1 - /3)/i/ = 0.87. Furthermore, it 
should be noted that the total variation with the lattice- 
size range is very small, ranging from l3/v — 0.56(1) to 
/3/v = 0.57(2), as we vary the size range from L = 8 — 44 
to L = 24 - 44. 

Summarizing, our findings in this Section on the 
emerging, under bond randomness, second-order phase 
transition of the 3d random-bond BC model are the fol- 



FIG. 12: FSS behavior of the order parameter at the esti- 
mated critical temperature. In the panel we show a simple 
power-law estimation of the exponent ratio /i/v. 



lowing: (i) the proposed critical exponents provide a sta- 
ble finite-size behavior, strongly supporting hyperscal- 
ing and (ii) the proposed value of the critical exponent 
'y/v — 1.864(12) characterizes, in a very clear way, the ex- 
pected distinctive strong-coupling fixed point, describing, 
according to the renormalization-group calculations, the 
emerging from the first-order regime second-order phase 
transition. 



V. SUMMARY AND CONCLUSIONS 

It is instructive at this point to attempt an overview 
on the effects of disorder for 3d ex-second- and ex-first- 
order phase transitions and compare our results with pre- 
vious relevant research. Let us restrict ourselves to a 
presentation focused mainly on some of the papers dis- 
cussed already in the text. These are the cases of the 
3d RIM [13, [13, [HI and those of the ex- weak first-order 
phase transition of the 3d site-diluted q = 3 Potts model 
studied by Ballesteros et al. [32] , and the case of the ex- 
strong first-order transition of the 3d bond-diluted q — 4 
Potts model studied by Chatelain et al. [SSJ . In table U 
we display the critical exponents obtained in these pa- 
pers together with the two concrete cases of the random- 
bond 3d BC Model. From this table one can see that, 
for the two cases studied here, the hyperscaling relation 
{2/3/v)+'j/v = dis well satisfied, and the best case is that 
of the ex- first-order transition that was also found to have 
a very robust FSS behavior. Furthermore, a straightfor- 
ward comparison shows that, our ex-second-order case 
{A — 1; r — 1/3) has a very similar behavior with that 
of Berche et al. [22] on the 3d bond-diluted Ising model 
with magnetic bond concentration p = 0.7. This is quite 
remarkable and, we may also point out that, our results 
are limited to lattice sizes up to L = 44, whereas sizes 
up to L = 96 have been simulated in that paper. Our 



12 



TABLE I: Summary of critical exponents for the 3d pure and disordered Ising (IM), g-states Potts (PM), and Blume-Capel 
(BCM) models, as obtained in Refs. [VA \2vl. |22. l32. 133. l64l. l65l| and the present paper. 



Model 




1/u 


V 


ll^ 


7) = 2--ilv 


filv 


{2filv) + 7/^ 


Ex- second- order phase transition 


Pure IM [64] 




1.586(3) 


0.6304(13) 


1.966(3) 


0.034(3) 


0.517(3) 


3.00(9) 


Pure IM [65] 






0.63002(10) 




0.03627(10) 






Site-diluted IM(p = 0.8)" [17] 




1.463(11) 


0.6837(53) 


1.963(5) 


0.037(5) 






Site-diluted IM(p = 0.8)" [20] 






0.683(2) 




0.036(1) 






Bond-diluted IM (p = 0.7)" [22] 




1.520(20) 


0.660(10) 


1.965(10) 


0.035(10) 


0.515(5) 


2.995(20) 


Bond-diluted IM (p = 0.55)" [22] 




1.460(20) 


0.685(10) 


1.977(10) 


0.023(10) 


0.513(5) 


3.003(20) 


Random-bond BCM (A = 1; r = 


1/3) 


1.504(19) 


0.665(17) 


1.964(4) 


0.036(4) 


0.510(10) 


2.984(24) 


Ex-first-order phase transition 


Bond-diluted g = 3 PM [32J 




1.449(10) 


0.690(5) 


1.922(4) 


0.078(4) 






Random-bond BCM (A = 2.9; r -- 


= 1/3) 


1.440(10) 


0.694(5) 


1.864(12) 


0.136(12) 


0.566(5) 


2.996(22) 


Bond-diluted g = 4 PM |33] 




1.330(25) 


0.752(14) 


1.500(14) 


0.500(14) 


0.645(24) 


2.790(62) 



^1 — p denotes the concentration of impurities. 



efforts indicate, in agreement with Berche tt al. [22|, that 
at these lattice sizes the system still crosses over to the 
universality class of the RIM, described by the second 
entry in table H] [13] • Finally, we have deliberately placed 
our second case of study of the random-bond BC Model 
(A = 2.9; r = 1/3) after the ex-weak- and before the ex- 
strong-first-order transitions of the 3d Potts model, since 
our results for the critical exponents appear to interpo- 
late between these two cases. 

In conclusion, the 3d random-bond BC model has been 
studied numerically in both its first- and second-order 
phase transition regimes by a comprehensive FSS anal- 
ysis. As expected on general universality arguments the 
3d random-bond BC model at the second-order regime 
(A = 1) was found to be fully compatible with the 3d 
Ising universality class. However, the case studied here 
[(A, r) = (1,1/3)] exhibits analogous crossover problems 
as the ones encountered in the case of the bond-diluted 
3d Ising model with magnetic bond concentration p = 0.7 
studied by Berche et al. [22| . For the case of ex- first-order 
regime at A = 2.9, we have shown that the disorder 
strength r = 1/3 was strong enough to convert, with- 
out doubt, the original first-order transition to a genuine 
second-order one. For this case [(A,r) = (2.9,1/3)] we 
have presented a detailed and convincing FSS scheme. 
The scenario adopted in this paper, and the proposed 
critical exponents obtained from a stable scaling behav- 
ior, are supporting strongly hyperscaling. In particu- 
lar the value of the critical exponent 7/1/ = 1.864(12) 
is robust and characterizes, in a very clear way, the ex- 
pected distinctive strong-coupling fixed point, describing, 
according to the rcnormalization-group calculations, the 
emerging from the first-order regime second-order pase 
transition. The present results point out, as in the case 
of the 2d random-bond BC model, the existence of a 
strong violation of universality principle of critical phe- 
nomena, since the two second-order transitions between 



the same ferromagnetic and paramagnetic phases have 
different sets of critical exponents. In the 2d random- 
bond BC model the difference in the exponents was re- 
vealed in the thermal exponent v and an extensive but 
weak universality was found corresponding to the fact 
that the two emerging transitions have the same mag- 
netic exponent ratios {P/v and ^/v) [il|. However, for 
the 3d random-bond BC model no such weak universality 
is supported. The observed strong violation of universal- 
ity is now revealed mainly, but not exclusively, in the fact 
that the corresponding emerging transitions have differ- 
ent magnetic exponent ratios, as seen from the values of 
•j/iy of table [D 

The proposal for the possible new universality class in 
the 3d random BC model in the ex-first-order regime 
is an interesting finding supported also by the early 
renormalization-group calculations [J, |lO[ , as mentioned 
already in the introduction. It is also a surprising result 
since the two transitions are between the same ferromag- 
netic and paramagnetic phases. However, having in mind 
that the 3d RIM suffers rather slowly decaying scaling 
corrections [13, 120, UM^ O'^s can never be confident for 
the asymptotic behavior at the present system sizes. In 
view of the above comments, and the earlier discussion of 
our fitting attempts for the susceptibility maxima, a con- 
fident resolution of this situation may require lattice sizes 
of the order of at least L = 240. Further investigations, 
for instance considering the behavior for different values 
of the crystal-field coupling A, suitably chosen in the ex- 
first-order regime, may be very useful to this direction. 
In a more advanced level one may try to obtain further 
convincing evidence by c omp aring multifractal spectra 
for correlations functions [gO, mM for the ex-second- and 
ex-first-order transitions. These lines of research demand 
further computational efforts and more sophisticated MC 
and FSS schemes, but as pointed out in subsection III Bl 
may yield more precise and additional useful information 
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concerning the properties of disorder averages \li 
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